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Influence of material stretchability on the equilibrium shape of 
a Mobius band 

David M. Kleiman," Denis F. Hinz,^ and Eliot Pried*^ 

We use a discrete, lattice-based model for two-dimensional materials to show that Mobius bands made with stretchable 
materials are less likely to crease or tear. This stems a delocalization of twisting strain that occurs if stretching is allowed. 
The associated low-energy configurations provide strategic target shapes for the guided assembly of nanometer and micron scale 
Mobius bands. To predict macroscopic band shapes for a given material, we establish a connection between stretchability and 
relevant continuum moduli, leading to insight regarding the practical feasibility of synthesizing Mobius bands from materials 
with continuum parameters that can be measured experimentally or estimated by upscale averaging. To take advantage of 
stretchability in the case of Mobius bands made of graphene, DNA, and other effectively unstretchable materials, we develop 
and explore a novel architecture that uses the Chinese finger trap as a fundamental building block and imparts notable 
stretchability to otherwise unstretchable materials. 


Recent technical advances have made it increasingly clear 
that the properties of a material are determined not only 
by its composition but also by geometrical and, impor¬ 
tantly, topological factors.® With this realization and 
breakthroughs in the ability to fabricate objects with 
molecular-scale precision, research into using the one¬ 
sided topology of the Mobius band in scientific appli¬ 
cations is burgeoning. In chemical topology, for ex¬ 
ample, mechanically interlinked molecules, or catanenes, 
have been created using Mobius molecules as interme¬ 
diaries, which may set the stage for synthesizing pro¬ 
grammable topological nanostructures.®^ A branch of 
algebraic topology known as homotopy figures promi¬ 
nently in the description of defects in ordered media.® 
In that context, every defect is characterized by a topo¬ 
logical charge, the conservation of which results in sta¬ 
bility. In particular, the stability of topological solitons 
and singular defects is connected to the topological sta¬ 
bility of Mobius bands.® An additional connection be¬ 
tween Mobius topology and defects is provided by mi¬ 
cron scale Mobius crystals, which were first created just 
over a decade ago by spooling niobium triselenide ribbons 
onto selenium droplets.® Such objects can be considered 
as global disclinations.® The ability of a recently synthe¬ 
sized expanded porphyrinoid to switch between Hiickel 
and Mobius top olog ies presents the possibility of novel 
memory devices.®^ Mobius topology is also exhibited by 
cyclotides, macrocylic plant proteins involved in plant de¬ 
fense.^ Due to their topologically derived structural sta¬ 
bility, these proteins have the potential to serve as drug 
scaffolds and pharmaceutical templates 

This article focuses on determining energetically pre¬ 
ferred equilibrium shapes of Mobius bands. Work on this 
problem was initiated by Sadowsky,®®^^^^^^^^^ who fo¬ 
cused exclusively on materials like paper which are easy 
to bend but essentially unstretchable and, thus, must 
adopt shapes that are very closely approximated by de¬ 
velopable surfaces. Aside from proving that it is pos¬ 
sible to construct a developable band from a rectangular 
strip of width sufficiently small relative to its length, Sad- 
owsky established an upper bound for bending energy of 
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a developable band and derived a dimensionally reduced 
expression for the bending energy of a band made from 
an infinitesimally thin rectangular strip. Wunderlichl^SIIIl 
later sharpened Sadowsky’s bound and generalized Sad- 
owsky’s bending energy to incorporate the effect of finite 
width. The problem of constructing developable equilib¬ 
rium configurations was first considered by Mahadevan 
and Kellerwhose numerically determined solutions led 
to a tighter upper bound on the bending energy but are 
inconsistent with results of Randrup and Rqgen,!^ who 
showed that the centerline of a Mobius band must have 
an odd number of switching points at which its curva¬ 
ture and torsion both vanish. Using Wunderlich’s energy, 
Starostin and van der Heijden®^^ computed equilibria 
that meet these requirements and also found evidence to 
suggest that Sadowsky’s energy is a singular limit that 
produces centerlines with discontinuous curvatures. 

We depart from established tradition and explore the 
influence of material stretchability on the shape of an 
equilibrated Mobius band. For bands with sufficiently 
large width-to-length ratios, Starostin and van der Hei- 
jden observed localized zones of concentrated bending- 
energy density. Their results point to the emergence of 
singularities indicative of the onset of failure. Even a 
slight degree of stretchability should alleviate such con¬ 
centrations. We seek to quantify both this effect and ac¬ 
companying shape variations that bands of a given aspect 
ratio exhibit with increasing stretchability. Additionally, 
a model that allows for small but discernible stretchabil¬ 
ity applies even to conventional paper, which stretches 
by a percent or two in the direction of loading without 
creasing or tearing.l^ 

Our approach utilizes a discrete, lattice-based model 
that incorporates stretching but can also accurately ap¬ 
proximate developable shapes for sufficiently small values 
of the stretchability. Regardless of the aspect ratio of 
the band, we find that the total energy decreases mono- 
tonicahy with stretchability. Bands made of stretchable 
materials are therefore easier to form than bands made of 
unstretchable materials. 

For further insight regarding how stretchability influ¬ 
ences equilibrium shape, we compute the mean and Gaus¬ 
sian curvatures for bands of various length-to-width as¬ 
pect ratios, the latter of which is nonzero only for stretch¬ 
able materials.® Consistent with the observation that un¬ 
stretchable materials must adopt developable shapes, we 
find that the Gaussian curvature plays a key role in bands 
comprised of such materials. Except for cases involv- 
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ing combinations of the smallest aspect ratio and the 
two largest values of stretchability investigated, the mean 
and Gaussian curvatures become more evenly distributed 
across the band as stretchability increases. Bending is 
concomitantly transferred to stretching, thereby eliminat¬ 
ing localized zones of concentrated bending-energy den¬ 
sity and the associated possibility of creasing or tearing. 

We also show that our model is energetically consistent 
with a simple continuum theory and derive relationships 
between our material parameters and those of the contin¬ 
uum theory. These relationships could provide a basis for 
future material design and applications. 

Finally, we describe a novel architecture consisting of 
a coiled Chinese finger trap. Our construction adopts a 
shape characteristic of those exhibited within the stretch¬ 
able regime of our lattice model. Comparing the particu¬ 
lar shape adopted by our construction with results from 
our model thus allows us to infer the effective stretcha¬ 
bility of the underlying architecture. 


1 Discrete, lattice-based model 

We first provide a brief description of our model. Details 
are provided in the Supplementary Information (SI). 

1.1 Kinematics 


1.3 Nondimensionalization 

The width re, length L, linear spring constant ki, and 
angular spring constant kg yield dimensionless parameters 

a = — and k = . ( 2 ) 

w w^ki 

While the (length-to-width) aspect ratio a embodies the 
geometry of the band, the stretchability k character¬ 
izes the ratio of bending resistance to stretching resis¬ 
tance. Smaller values of k describe materials like paper, 
graphene, and DNA which bend easily but are difficult 
to stretch. The limiting case k ^ 0 of vanishing stretch¬ 
ability embodies the idealized limit of a material which 
cannot be stretched and thus can adopt only developable 
shapes. 

Next, using A = Lw as a reference area, kiA as a ref¬ 
erence energy, and defining 

E , r(i) . ~ Te 

r(t) = -, and = — (3) 

kiA w w 

leads to a dimensionless version 

«' = «'« + '^b, (4) 

of the total energy 0 , with stretching and bending con¬ 
tributions 


The shape of an unstretchable Mobius band is uniquely 
determined (up to a rigid transformati on) by the cur¬ 
vature K and torsion r of its centerlineThe shapes 
adopted by such bands must be developable. That is, 
they must be ruled surfaces that can be continuously flat¬ 
tened onto planar regions without stretching/contracting 
or, equivalently, with the preservation of intrinsic lengths 
and angles. In contrast, a band made of a stretchable 
material can adopt a nondevelopable equilibrium config¬ 
uration, in which case its shape is uniquely determined 
(up to a rigid transformation) by its first and second fun¬ 
damental forms, as discussed in the SI. 


1.2 Linear and angular springs 

We approximate a rectangular strip of length L and width 
re by a lattice of equilateral triangles with N points uni¬ 
formly separated by a distance vq. To incorporate re¬ 
sistance to stretching, we connect each pair of lattice 
points by a linear spring with stiffness ki and equilibrium 
length Te- Further, to incorporate resistance to out-of¬ 
plane bending, we connect each triplet of lattice points 
with a torsional spring of stiffness ko > 0 and equilibrium 
angle Oe = tt. The total energy E of the band is then 
given by the sum 

E = ^ + T (1) 

ieSi ieSe 


where Si and So denote the sets of all linear and torsional 
springs, r{i) is the current length of the i-th linear spring, 
and 0{i) is the angle between triplets of points associated 
with the i-th torsional spring. 




E 

ieSi 


(r(i)-fe)^ 

2a 


and T 5 


E 

ieSo 


For brevity, the dimensionless quantities 
are hereafter referred to as energies. 


kjeji) - 0ef 
2a 

(5) 

T, T,, and T 5 


2 Simulation results 

Our simulations indicate that stretchable Mobius bands 
adopt three characteristic equilibrium shapes, depending 
on the combination of stretchability k and aspect ratio 
a. For small values of k, the bands take shapes that 
resemble those of paper models and previous simulation 
results. Increasing k leads to a loss of developability, and, 
for a certain combination of low aspect ratio and high 
stretchability, a self-intersecting achiral shape emerges. 

2.1 Equilibrium shapes 

Equilibrium band shapes obtained by minimizing the to¬ 
tal energy T are provided in Figure for representative 
combinations of k and a. For sufficiently small values 
of k and each value of a considered, equilibrium shapes 
qualitatively resemble those of model bands made from 
rectangular strips of paper and thus appear to be nearly 
developable. Regardless of the value of a, reducing k 
appears to have a negligible influence below k = 10 “^. 
However, the influence of k becomes progressively more 
evident above k = 10 “^ and is increasingly obvious for 
smaller values of a. The centerlines of bands appear to 
be more circular and less out of plane for larger values of 
/c, an impression that is confirmed by plots for bands of 
aspect ratio a = 27r that appear in Figure 
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Figure 1: Mobius bands adopt characteristic equilibrium shapes depending on stretchability, fc, and aspect ratio, a: Eq uilibr ium shapes 
for different values of a and k. Bands are rotated into their main axes. For views from other perspectives, see SI Figure [SI.5| 


If k is sufficiently small, the tangent vector to the cen¬ 
terline of each equilibrated band exhibits an odd number 
of switching points, at which its curvature, k, and torsion, 
r, vanish simultaneously. While k has two peaks and one 
zero, r has two peaks and three zeroes (Figure]^. These 
findings are co nsistent with previous analytical and nu¬ 
merical result Moreover, for smaller values of a, 
the peak values of k and r are higher. In keeping with 
findings of Mah adeva n and KelleiSSl and Starostin and 
van der Heijden,II2l2] centerlines of these bands are 
more out-of-plane than those of bands with larger aspect 
ratios (Figure®. See the SI for comprehensive conver¬ 
gence and valiaation studies in the unstretchable limit 



Figure 2: Centerlines of the equilibrium shape of bands with aspect 
ratio a = 7T made from materials of various stretchabilities k. 



Figure 3: In the unstretchable limit, the model recovers analytical 
predictions of Randrup and Rpgenl^ and t he cha racteristic shape 
obtained by Starostin and van der Heijden. h^ l i^ l Curvature k and 
torsion r of the centerline of Mobius bands for stretchability k = 
10 and various values of the aspect ratio a versus the arclength 
along the centerline in terms of the number Ni of points along the 
centerline. 


k^O. 

For a = IT^ bands made of materials with stretchabil¬ 
ities k = 10“^ and k = exhibit self-intersecting 

achiral equilibrium shapes. This degeneration, which is 
possible only if unpenalized self-intersections are allowed, 
resembles the collapse observed in Mobius soap films with 
small throat distances.!^ However, the constraint of lat¬ 
tice connectivity inherent to our model delivers shapes 
different from those of collapsed Mobius soap films. Ad- 
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Figure 4: The shape of the centerline of an effectively unstretchable 
Mobius band depends on its aspect ratio a: Centerlines of equi¬ 
librium shapes of Mobius bands for stretchability k = and 

representative values of a. Bands are rotated into their main axes 
as in Figure^ 


ditionally, the choice k = 10 yields not a band but a 
flat annulus. 


gineering. The synthesis of Mobius bands is challenging 
mainly due to their significantly higher energy states rel¬ 
ative to untwisted rings. However, allowing for stretching 
diminishes this energy difference and could therefore re¬ 
duce the difficulty of fabrication strategies that rely on 
bending rectangular strips. 

Since T decreases as k increases for each value of a 
considered, the loss of chirality exhibited by bands of as¬ 
pect ratio a = TT made from materials with stretchabil- 
ities k = 10“^ and k = 10^ indicates that chirality is 
is accompanied by increased stored energy. To fabricate 
chiral objects, a system must therefore be sufficiently and 
properly constrained. For example, anisotropic particles 
can be used to guide the assembly of chiral objects.l^ 
However, we speculate that, for k sufficiently large, non- 
collapsed bands are at best metastable, since exhibits 
oscillations at k = 10^ observed for small values of a. 
If a material is too stretchable, then small aspect ratio 
Mobius bands are expected to be unstable. This provides 
a threshold for the stretchability below which the synthe¬ 
sis or guided assembly of small aspect ratio Mobius bands 
should become feasible. 


2.2 Energy measures 

Plots of the stretching energy the bending energy 
and the total energy T, all normalized by k to en¬ 
compass changes in overall stiffness, are provided in Fig¬ 
ure for representative combinations of k and a. For 
each choice of a, first increases as k increases but is 
dominated by the order-of-magnitude larger which 
decreases. Hence, the sum T = Tg + decreases mono- 
tonically as k increases and the minimum value of T for 
any band is attained at the largest value of k considered. 
From an energetic perspective, it therefore seems reason¬ 
able to infer that bands made of stretchable materials 
are significantly less costly to make than bands made of 
unstretchable materials. 

The influence of a is most evident for k less than a 
critical value kc ^ 10 “^. For each k less than kc^ T 
decreases monotonically with a. This effect is in keep¬ 
ing with everyday experience: the effort needed to twist 
a rectangular sheet of an essentially unstretchable mate¬ 
rial like paper into a Mobius band increases notably as 
its length-to-width aspect ratio deceases. Moreover, per¬ 
forming the same task with a stretchable material, such 
as a thin sheet of rubber, is substantially easier. Inter¬ 
estingly, however, for k > kc and a sufficiently small, T 
appears to increase, albeit weakly, with increasing a. 

With reference to the previously noted observation that 
the centerlines of bands become more circular and less 
out-of-plane as k increases, the observed dependence of T 
on k indicates that more planar, less bent configurations 
are energetically favored above a certain critical stretcha¬ 
bility. For the particular choice k = 10 “^, is negligibly 
small in comparison to and, thus, T ^ The choice 
k = 10“^ therefore suffices to capture the strictly un¬ 
stretchable limit k = 0 . The influence of k is otherwise 
most significant near that limit. This provides potentially 
important insight regarding the guided assembly or syn¬ 
thesis of Mobius bands, as it indicates that small increases 
in stretchability yield significant gains in the ease of en¬ 


2.3 Curvature and dilatation 

In the SI, we show that the coarse-grained limit of our 
model corresponds to a continuum model for a surface S 



Figure 5: Stretching energy Tg, bending energy and total en¬ 
ergy T for representative combinations of stretchability, k, and as¬ 
pect ratio, a. These quantities are normalized by k to encompass 
changes in overall stiffness. 
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Figure 6: The mean and Gaussian curvatures are both more evenly distributed for increasing stretchability k: (a) Square of the mean 
curvature if, (b) magnitude \K\ of the Gaussian curvature iF, and (c) square of the strain e of Mobius bands made of stretchable 
materials for different values of the aspect ratio a and stretchability k. The aspect ratio of each contourplot is equal to the aspect ratio 
of the band it represents. Images of the bands corresponding to the contour plots appear in Figur e [T] N ote that different scales are 
adapted individually to capture the entire range of relevant values. Additional cases appear in Figure [SI. 6 1 of the SI. 


that resists stretching and bending. Whereas stretching is 
characterized by an area modulus fia related to the linear 
spring stiffness ki by 


With ^ and 0, the total energy £ of S takes the form 

S = Vs [ [kie^ + - K)] da, (8) 

Js 


l^a = 


Vski 

2 


(6) 


bending is characterized by splay and saddle-splay moduli 
IV and jv related to the torsional spring stiffness kg by 


Sy^kg ^ _ 2/i ^ 

11 =—-— and ii = -^ = V 3 kg. ( 7 ) 


where da denotes the area element on S and where 
i 7 , and K denote the (two-dimensional) dilatation, mean 
curvature, and Gaussian curvature of S. The expression 
® for the area modulus fia is consistent with a result 
Stained by Seung and Nelson.l^ The bending contribu¬ 
tion to 0 arising from Q, which is identical to that of a 
Kichhoff plate with bending modulus 3 y^kg /2 and Pois¬ 
son’s ratio 1 / 3 , coincides with an expression derived by 
Merchant and Keller 
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Figure 7: Both mean and Gaussian curvature are concentrated along the intersection of self-intersecting achiral bands: (a) Collapsed 
bands colored by stretchability as in Figure(b) Square of the mean curvature 77, (c) magnitude |7F| of the Gaussian curvature 7F, 
and (d) square of the strain e of collapsed Mobius bands for different stretchabilities k. The aspect ratio of each contour plot is equal 
to the aspect ratio of the band it represents. Note that different scales are adapted individually to capture the entire range of relevant 
values. 


This connection to continuum theory suggests that the 
pointwise distributions of suitable discrete approxima¬ 
tions to and may provide further insight regard¬ 

ing how the stretchability k and aspect ratio a influence 
the shape of a Mobius band. Plots of these distributions 
are provided in Figure for representative values of k 
and a. In the effectively unstretchable case k = 10 “^, 
each equilibrium shape exhibits a nearly flat triangular 
region bounded at its vertices by zones in which takes 
large values and is surrounded by a nearly flat trapezoidal 
region. These ch aracterisics are consistent with previous 
findingsas is the increase in size of the zones in 
which exhibits large values with increasing a. How¬ 
ever, as k increases, becomes increasingly more evenly 
distributed over the band. Consistent with the observa¬ 
tion that localized regions of concentrated bending-e nergy 
density indicate where creasing or tearing may occur ,^^11 
our findings show that for bands of aspect ratio increasing 
the stretchability alleviates such concentrations. 

All equilibrated bands exhibit non-vanishing Gaussian 
curvature K. However, for k = 10 “^, K is very close to 
zero almost everywhere. This confirms that the choice 
k = 10“^ yields nearly developable shapes and thus pro¬ 
vides a good approximation of the unstretchable limit, 
even though some degree of stretching is allowed for any 
value of k > 0 . For all larger values of k considered, 
and \K\ are maximized at the same points. The maxi¬ 
mum values of I if I become more prominent as a decreases. 
At the vertices of the previously discussed flat, triangu¬ 
lar regions, where the contribution to the bending energy 
from is largest, the contribution to the bending energy 


from \K\ also attains its largest values. Thus, even in the 
approximately unstretchable limit, it is energetically fa¬ 
vorable to locally transfer energy associated with to 
energy associated with K. Inspecting the local dilatation 

shows that bending is locally transferred to stretching 
unless K = 0 , in which case the shape of the band is gen¬ 
uinely developable. Since most materials manifest some 
in-plane elasticity, our findings also indicate that an ac¬ 
tual band should possess a region of local stretching. This 
also suggests that loci of large are likely to coincide 
with zones in which the bending contribution to the en¬ 
ergy is concentrated. These zones are also most likely to 
be loci for inelastic deformation or failure. Conversely, 
allowing for local stretching may reduce the bending con¬ 
tribution to the energy density. 

Both and \K\ become more evenly distributed over 
the band with increasing /c, resulting in lower curvature 
gradients. Further, both \K\ and are transferred from 
the edge toward the centerline of the band, reducing the 
magnitude of the gradient of the continuum bending- 
energy density and decreasing the likelihood of creasing or 
tearing accordingly. Nonzero values of |i^| indicate that 
stretching contributes significantly to the overall shape 
of the band, as zones in which \K\ is largest coincide 
with zones of large , to which the continuum stretching- 
energy density (see the SI) is directly related. This result 
confirms the previous observation that increasing k causes 
the centerlines of bands to become more circular and less 
out-of-plane, thereby reducing bending. A zone of low 
bending energy is maintained as k increases. The magni¬ 
tude \K\is largest in zones in which takes large values. 
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Figure 8: Lines of curvature are symmetric and trace out a nearly 
regular grid for bands made of stretchable materials, indicating re¬ 
duced concentrations of bending-energy density: Lines of curvature 
for three representative characteristic shapes are shown on a rectan¬ 
gular strip with the associated equilibrium shape provided above, 
(a) Developable shapes have pronounced lines of curvature along 
the length of the band with nearly straight lines of curvature along 
their widths, (b) For bands made of stretchable materials, shapes 
with circular centerlines also have nearly straight lines of curvature, 
however any such band curves more gradually near its axis of sym¬ 
metry. (c) Collapsed shapes have nearly straight lines of curvature 
with a sharp curve near the center of the intersection. 


resulting in large values of the continuum bending-energy 
density. However, this effect is mitigated by reduced gra¬ 
dients due to stretching. 

The transition from bending to stretching occurs 
smoothly except for combinations of relatively large 
stretchability k and very low aspect ratio a leading to col¬ 
lapsed, achiral equilibrium shapes. Plots of |i^|, and 
are provided in Figure for three illustrative choices 
of k and a = it. For k = 10“^ and k = and 

\K\ appear to diverge at a single point along the inter¬ 
section and decay quickly with increasing distance from 
that point. The value of is also greatest at the point 
in question but does not decay as rapidly. However, for 
the choice k = which results in collapse to an an¬ 
nulus, the values of \K\^ and are largest along the 

centerline of the band, decaying toward the edge of the 
band. 


2.4 Lines of curvature 

Differences in the distributions of curvature and stretch 
exhibited by the three characteristic shapes of Mobius 
bands are also evident from the lines of curvature, which 
are curves with tangent vectors that align with a prin¬ 
cipal direction at every point. Plots of these curves for 
representative combinations of stretchability, k, and the 
aspect ratio, a, are provided in Figure 

The least regular grid occurs for k = 10“^. The im¬ 
provement in regularity that occurs as k increases indi¬ 
cates that smaller gradients of the curvature correspond 
to smaller strain gradients. In agreement with the previ¬ 
ous discussion, we thus infer that bands made of stretch¬ 
able materials have reduced concentrations of bending. 


The lines of curvature of a collapsed band form a nearly 
regular grid but curve sharply near the intersection of the 
band (which coincides with the previously mentioned axis 
of symmetry). In the absence of collapse, for bands made 
of stretchable materials we otherwise observe similar be¬ 
havior. However, lines of curvature bordering the axis of 
symmetry of such a band are smoother. 


3 Applications 


3.1 Graphene 

Based on available data for the (two-dimensional) 
Young’s modulus Y and (two-dimensional) Poisson’s ratio 
u of graphene, its area modu lus /Tq = Y/2(l+z^) should be 
no less than 130 N/Using this lower bound in ^ 
shows that, to model graphene, the linear spring stiffness 
must satisfy 

ki > 150 N/m. (9) 

In contrast to the classical Kirchhoff theory, the bend¬ 
ing moduli of graphene are independent of Y and uJ^Foi 
the reported value 22 • 10“^^ Nm of the bending modulus 
of graphene,l^^l^^the torsional spring stiffness determined 
by 01 is 

4 = 85 • 10“^^ Nm. (10) 


Granted that ki and ko satisfy 0 and (10), we may 
use the definition ( 1^2 of the stretch^ility k to determine 


upper bounds on the width w along with associated lower 
bounds on the length, L, of a graphene band of given 
aspect ratio, a > tt, needed to give rise to a targeted 
value of k. To achieve k < 10“^, for instance, it would be 
necessary to take re < 2.4 nm (corresponding to a width of 
approximately 10 unit cells) and L > 2.4a nm. Reducing 
the stretchability by another order of magnitude would 
require reducing the width to a few unit cells. Except for 
exceptionally large values of a, the equilibrium shape of a 
Mobius nanoband made from graphene should therefore 
be nearly developable. 


3.2 A stretchable material architecture 

As a consequence of the structurally induced effective 
stretchability of a discrete lattice, there is a possibil¬ 
ity of using a nanomaterial with novel architecture to 
fabricate Mobius bands without regions of concentrated 
bending-energy density. As a particular example, effec¬ 
tively unstretchable graphene nanoribbons (GNRs) could 
be formed into a stretchable biaxial braid with the struc¬ 
ture of a (Chinese) finger trap and then assembled side- 
by-side to form a ribbon. 

Such ribbons could then be twisted into a Mobius 
band. Since our calculations as well as the work of Qi 
et show that the behavior of strained graphene at 

the nanoscale is qualitatively similar to that of strained 
paper at the macroscale, it is instructive to realize this 
construction with a paper model. The macroscale Mobius 
band shown in panels (c), (d), and (f) of Figure [^adopts 
a shape close to that characteristic of bands with stretch¬ 
ability k = 10“^ and aspect ratios satisfying a > tt. On 
this basis, we estimate that the stretchability of our pro¬ 
posed architecture is roughly equal to k ^ 10“^. More 
generally, this example consequently suggests that our 
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Figure 9: A Mobius band formed from finger trap structures adopts a characteristic stretchable shape: (a) One strand of the finger trap 
structure, (b) The three strands glued together, (c-g) Stretchable finger trap Mobius band from different perspectives compared to 
simulation results for a = 27 t and k = 10“^. 


model can be used to infer the effective stretchability of 
a material from the shape it adopts when twisted into a 
Mobius band. 

DNA origami techniques pioneered by Han and cowork- 
ersEl might be used to fabricate nanometer or micron scale 
versions of the macroscopic finger trap Mobius bands con¬ 
sidered here. 

4 Conclusions and discussion 

Using a discrete, lattice-based model, we confirm analyti¬ 
cal results for the distributions of curvature and torsion on 
the centerlines of Mobius bands^^ along with results qual¬ 
itatively similar to those obtained for con tinuu m mod- 
eis [MEl analogy to the effect of writhe,l^^^^ we find 
that twisting strain becomes delocalized if stretching is 
allowed. To predict macroscopic band shapes for a given 
material, we establish a connection between stretchabil¬ 
ity and relevant continuum moduli. This affords insight 
regarding the practical feasibility of synthesizing Mobius 
bands from materials with continuum parameters that 
can be measured experimentally or estimated by upscale 
averaging. 

We find that Mobius bands adopt three characteristic 
shapes depending on stretchability, k, and aspect ratio, 
a. These provide reasonable target shapes and strategies 
for the guided assembly and synthesis of Mobius bands. 
In particular, we find that highly stretchable bands have 
lower energies, indicating that it would be easier to make 
Mobius bands from more stretchable materials. Hence, 
Mobius molecules might well correspond to these stretch¬ 
able shapes. We also find that for bands with k < 10“^, 
the smaller a becomes, the higher the energy, so that 
wider bands should be more difficult to make than nar¬ 
rower ones. Conversely, we find that for k > 10“^ bands 
with lower a have lower energies, but for a = tt and either 
k = 10“^ or k = they collapse into self-intersecting 

achiral shapes, which are degenerate Mobius bands. This 
suggests that highly stretchable bands with small aspect 
ratios might be unstable and thus impossible to make. 
For 10“^ ^ k < 10“^, bands with smaller values of a 


should, however, be easier to make than those with larger 
ones. 

Thus, k and a influence not only the equilibrium shape 
of a Mobius band but also its ease of guided assem¬ 
bly. Due to the energetically advantageous nature of 
bands made of stretchable materials, the ability to design 
and fabricate Mobius bands with nondevelopable shapes 
might be valuable. These findings could be helpful in a 
variety of applications and also serve as an archetype for 
other twisted topologies. 

In particular, our results indicate that bands made of 
stretchable materials are distinguished from bands with 
developable shapes by lower curvature gradients. Conse¬ 
quently, the stretchability k is positively related to the 
extent to which the curvature of a band is homogeneous. 
This connection between stretchability, curvature gradi¬ 
ents, and shape has implications for the electronic ground 
state of quantum Mobius bands .1^ In particular, the 
ground and excited states of a Mobius band determined 
by minimizing Wunderlich’s functional over a restricted 
family of centerlines were found to exhibit wave functions 
that are significantly altered by curvature effects.l^ Sim¬ 
ilarly, curvature has been found to influence electron lo¬ 
calization, where deep potential wells arise in connection 
with singularities of the bending-energy density.l^ Since 
stretchability and curvature are inextricably linked,!^ our 
results suggest that the electron localization of a quantum 
Mobius band is directly related to its stretchability and 
aspect ratio. For bands made of stretchable materials, 
the combined influences of nontrivial Gaussian curvature 
and mean curvature may yield unprecedented effects. 

Mobius bands comprised of graphene with architectures 
that impart stretchability might allow for advances based 
on Mobius topology. Importantly, the use of graphene 
Mobius bands for topological insulators has already been 
invest igated.l^ Here, we presented a simple but realistic 
finger trap paper model as macroscopic analog for stretch¬ 
able graphene Mobius structures that allow determina¬ 
tion of the underlying stretchability. Such analog models 
could thus be used as a straightforward means to test new 
architectures and gain further insight into their proper- 
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ties. Further, couiparing the shape adopted by our finger 
trap Mobius bands with numerical results generated by 
our model yields a quantitative estimate for the effective 
stretchability of the novel woven construction proposed 
here. 

To suppress the fragility and fast decoherence of quan¬ 
tum states and yield topological stability, a notion resem¬ 
bling that underlying our finger trap model has been pro¬ 
posed for the formation of a topological qubit by braiding 
Marjoranas.S^®^ Also in the quantum realm and taking 
advantage of Mobius topology, a method to enhance spi¬ 
ral intramolecular charge transfer in Mobius cyclacene for 
use in novel optical and photoelectric devices has been 
proposed.l^ 

Given the ever-increasing spectrum of applications, the 
treatment of stretchability and aspect ratio and their im¬ 
plications for the Mobius band presented here may help 
guide further developments as well as to increase the the¬ 
oretical understanding of these fascinating objects. 


5 Methods 


5.1 Energy minimization 

To minimize the total energy T appearing in Q , we use a 
conjugate-gradient method within the molecular dynam¬ 
ics code LAMMPS (Large-scale Atomic/Molecular Mas¬ 
sively Parallel Simulator) .SSI For all values of the aspect 
ratio a considered, trial configurations for computing en¬ 
ergy minima are provided by the piecewise isor netri c con¬ 
struction (Figure 10) described by Sadowsky.SESl Since 
each trial shape is onvelopable, minimizing the energy of 



Figure 10: Schematic of a trial configuration provided by Sad- 
owsky’^lE^ developable Mobius band and used for computing the 
energetically preferred shape of a band of aspect ratio a = Gtt dis¬ 
cretized by N = 1170 points separated by uniform dimensionless 
distance rojL = \/3/727r. 


a band made of a material with stretchability k very close 
to zero constrains stretching and thus maintains approx¬ 
imate developability. With the linear springs already at 
their equilibrium distance, this minimization amounts to 
reducing the curvature of the band without altering the 
distances between points. Conversely, for non-negligible 
values of k, the spacing between points is unconstrained. 
The model and energy minimization strategy were sub¬ 
jected to extensive validation procedures, as described in 
the SI. 

5.2 Curvature approximations 

At each discrete point on the surface of the band, we com¬ 
pute discrete versions of the mean curvature iL, the Gaus¬ 
sian curvature K, and the principal curvatures ki and k ,2 
from discrete versions of the first and second fundamental 
forms. To calculate the latter quantities, we proceed by 
analogy to the continuous case^ where tangent vectors 
are approximated by vectors between neighboring points 
on the surface. 

5.3 Lines of curvature 

We determine the vector field of principal directions and 
approximate the lines of curvature by computing the rel¬ 
evant streamlines. All computed lines of curvature are 
symmetric about an axis perpendicular to the length of 
the band. 


5.4 Finger trap models 

Paper models for the finger trap Mobius band were made 
using standard office paper (80 g/m^) and liquid craft 
glue. 
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6 Discrete approximations of various curvatures 


We approximate a rectangular strip of length 
uniformly separated by a distance vq. Th e pie ce 
bend the strip into a Mobius band (Fignre |SI.l| ), 

cn plam 


and width re by a lattice of equilateral triangles with N points 
iiecewise developable mapping described by Sadowskyl^ ig then used to 
which we then employ as the trial configuration for each numerically 
dete rmined band with aspect ratio a. Each plane strip is discretized with a lattice of equilateral triangles (SI Figure 
SL2). The variable number of points used to discretize the lattice is denoted by N. Further, the variable number of 


points along the centerline is denoted by Ni. 

For a discretized curve with point spacing rp, the discrete curvature K{i) at the i-th point along the curve is given 

«(i) = ’'-^W')-*(‘ + i)) . (SI.1) 


^0 


where Z(t(i),t(i + 1)) denotes the angle between discrete tangent vectors at the point i —i.e., the angle at point i 
between points i — and i + 1 (see, for example, Belyaevl^. Similarly the discrete torsion r(i) at the i-th point 
along the curve is given by 


^ Z(b(i),b(z + 1)) 
ro 


(SL2) 


where Z(b(i), b(i + 1)) denotes the angle between the discrete binormals to the curve at points i and i + 1. 

To calculate the first and second fundamental forms of surface represented by a discrete triangular lattice, we 
proceed by analogy to the continuous case (see, for example, do CarmcP), but we use vectors between neighboring 
points on the surface as tangent coordinate vectors. At each point i of the lattice, we determine the discrete mean 
curvature i^(i), the discrete Gaussian curvature K{i), and the discrete principal curvatures Ei{i) and /^ 2 (i) from the 
first and second fundamental forms. 


y 



Fig. SI.l: Sadowsky’s parametrization of a Mobius band with a developable shape in a Cartesian coordinate system colored 
by the z-coordinate of each point. 


1 











Fig. SL2: A lattice of equilateral triangles with the three directions of angular and linear springs illustrated. 


Let a be the angle between the tangent vector to a curve and the principal direction associated with ki. The 
discrete normal curvature at the i-th point along a curve can then be determined pointwise from the continuous 
relation 

Kn = 1^1 cos^ a hv2 sin^ (SI.3) 


Further, having determined K{i) and /^n(^) from (SI.l) and (SI.3), 
curvature hig{i) can be determined from the continuous relation 


the magnitude \K,g{i)\ of the discrete geodesic 


\Kg\ = kI. 


(SL4) 
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Fig. SL3: The model shows convergence for increasing number of points along the midline: The curvature n (a) and torsion 
T (b) of the centerline for approximately unstretchable bands (k — 10“®) approach a limiting value for a large number of 
points along the centerline 173 < Ni. Notice that, as Ni is further increased beyond Ni = 346, slight oscillations in n and r are 
observed, indicating an ideal range of 173 < Ni < 346. The total energy T (c), stretching energy Tg (d), and Tj, (e) converge 
to a limiting shape corresponding to the continuum limit. Results are shown for a = 47r. 
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Fig. SL4: The centerline of a approximately unstretchable Mobius band depends on its aspect ratio: Centerlines of equilibrium 
shapes of general Mobius bands for different values of the aspect ratio a and stretchability k. The band is rotated into its 
main axes. 
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Fig. SI.5: Generalized Mobius bands adopt characteristic equilibrium shapes depending on the in-plane stretchability of the 
surface and the aspect ratio: Equilibrium shapes of generalized Mobius bands for different aspect ratio a and stretchability k. 
The band is rotated into its main axes. 



(a)fc=10° fc = 10"^ fc = 10"^ fc = 10"® fc=10“^ fc = 10“® fc = 10"® 



(b) fc = 10° k = 10"^ k = 10"^ k = 10"° k = 10“^ k = 10“° k = 10"® 



(c) fc = 10° k = 10“^ k = 10“^ k = 10“® k = 10“^ k = 10“® k = 10“® 



Fig. SI.6: Both mean and Gaussian curvature are more evenly distributed for increasing stretchability k: (a) Square of 
the mean curvature H, (b) magnitude |/vr| of the Gaussian curvature, and (c) square of the strain e of generalized Mobius 
bands for different values of the aspect ratio a and stretchability k. The aspect ratio of each contour plot is equal to the aspect 
ratio of the band it represents. Note the different scales adapted individually to capture the entire range of values. 
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7 Model validation 


To validate our model, we first investigate the influence of the discretization on the equilibrium shape and correspond¬ 
ing energy obtained after minimization. In view of the preceeding discussion, we expect that the model approaches 
a continuum limit as the number N of mesh points is increased. We study the impact of mesh size. To that end, we 
consider the convergence of the discrete approximations of the curvature k and the torsion r of the centerlines for 
bands made from effectively unstretchable materials, whereas for Mobius bands made from stretchable materials we 
study the convergence of T 5 , and T. 

Since k and r respectively depend on the directions of t and n, given a point spacing vq their convergence does 
not depend on a length scale. Instead it depends only on the number of points along the centerline. If we denote 
the arclength of the centerline by 5 , the arclength at the i-th point along the centerline by and the number of 
points along the centerline of the band by Ni, then it is most suitable to compare the discrete approximations of 
the curvature and torsion of the cent erline for different values of Ni. These approximations are observed to converge 
to limits as Ni is increased (Figure |SI.3[a-b)). Our results indicate that the shapes of bands made of effectively 
unstretchable materials remain unchanged on increasing Ni from 173 to 347, demonstrating that a choice of Ni in 
the range 173 < Ni < 34 7 prod uces converged results and that further increasing Ni has negligible impact on the 
shape of the band (Figure [SL3| ^a-b) ). 

For bands made from stretchable materials, we comp are T , and T 5 for increasing Ni and different values of k. 
All energies show convergence to limiting values (Figu re SI.3T-e)), which is consistent with the previously discussed 
convergence of the curvature and torsion (Figure [SI. 3 ^a-b)) in the case of bands made from effectively unstretchable 


materials. Note that the value of k for which energy due to stretching becomes negligible and the bands approach 
the unstretchable limit is not immediately evident. However, for k = 10“^, Tg is approximately 1% of T, indicating 
that it is reasonable to conclude that the choice k = 10 “^ provides a very close approximation to unstretchability. 


8 Connection between lattice and continuum models 


8.1 Bending energy 

Sadowskywork on Mobius bands made from effectively unstretchable materials hinges on a continuum descrip¬ 
tion in which a strip is treated as a surface S endowed with bending-energy density proportional to the square of its 
mean curvature H. By the Theorema Egregium of Gauss,bending a surface without stretching or shrinking leaves 
its Gaussian curvature K unchanged. Simultaneous bending and stretching of a surface will, however, generally alter 
its Gaussian curvature. To model stretchable bands, it is therefore natural to incorporate energetic dependence on 
Gaussian curvature. In keeping with a quadratic dependence on mean curvature, this amounts to adding a term 
linear in the Gaussian curvature and, granted that the band is free of spontaneous mean curvature, leads to the 
Canham-Helfrich bending-energy functional 


-^bend 


f 

Js 


^+T1K) dA, 


(SL5) 


involving the bending-energy density 


t^bend — T jlK 


(SL 6 ) 

where /i > 0 i s the splay modulus and /i < 0 is the saddle-splay modulus. Importantly, sinc e K = 0 for a developable 
surface, (SI.5) reduces in the unstretchable limit to the functional consid ered by Sadowsky^^^^. 

Although connections between the Canhar^Helfrich functional (SI.5) aiM triangular lattice parameters have al¬ 
ready been established by Seung and Nelson^^ and Schmidt and Fraternalil^, these rely on potentials involving the 
normals of neighboring triangular cells. Defining the normals to the cells of a triangulation of a nonorientable surface 
like a Mo bius b and is problematic. Hence, we use other elementary geometric considerations to establish a connection 
between (SI.5) and the bending-energy functional 


E, 




(SI.7) 


ieSe 


ieSe 


of the discrete lattice, where So is the set of all angular springs. 

To determine the effective spla y and saddle-splay moduli fi and ft corresponding to the underlying triangular 
lattice, we first use the definition (SI.l) of the discrete curvature hi{i) in the expression for the bending energy E^. 
Taking into account that the lattice is periodic lattice, this yields 


N 


Eb = 


(SL 8 ) 


i=l j=l 


where vq is the distance between lattice points, N is the number of points on the surface, and ko > 0 is the lattice 
modulus associated with bending {Eo in (SI. 8 ) carries the same dimensions as F^bend in (SI.5)), is the curvature 
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of the j-direction lattice line at point i for j = 1, 2, 3. Further, letting the angle betwee n the 1-dire ction lattice line 
at point i and th e first principal direction at point i be (a(i), we can use the expressions (SL3) and (SL4) for and 
n to write (SL 8 ) as 


N / 3 

Eb = \k0rl ^ + cos^{a{i) + b + cos^(a(z) + ^)) 


i=l \i=l 


3 / ' ^ ^ ^ ' 3 

+(K2(0)^(sin^(o;(0) + sin^(a(z) + b + sin^(a(z) + ^)) ) , (SI.9) 


where the angle between lattice lines is required to remain equal to 7 r/ 3 . To enforce this constraint on the angle, we 
assume the distance between points remains constant. This assumption yields reasonable approximations for small 
deformations. Further, it implies that lattice lines are geodesics—i.e., Kg = ^ along lattice lines. We next use the 
identities 

TT Q 

cos^ 'd + cos^('d + —) + cos^('d + -—) = -, 


TT Q 

sin'^^ + sin'^(^ + -) +sin'^(^ + —) = -, 
3 3 o 

cos^ + cos^(t? + ^) + cos^(i9 + 

o o Z 

sin^ + sin2(i9 + b + sin2(i9 + ^) = 5 

3 o Z 


(SI.IO) 


to reduce (SI.9) to 


Eb = 




2 ^ V 2 

i=l 


(SI.ll) 


where H{i) and K{i) are the respective discrete mean and Gaussian curvatures at the lattice point i. Finally, a 
Riemann sum argument yields 


Eb « Vske [ - K) 

Js 

for which the corresponding bending-energy density is 

^hend = VSkeiSH^ - K). 


Comparison with (SI.5) yields 


3\^ko ^ - /^7 

/i = —-— and fi = —voko. 


(SI. 12 ) 


(SI.13) 


(SI.14) 


To validate the foregoing calculations, we pres ent th e example of a plane strip discretized with a lattice mapped 
isometrically to a helically wrapped strip (Figure SL 2 ). Relative to a rectangular Cartesian basis { 61 , 62 , 03 }, 

Let { 01 , 62 , 63 } denote a rectangular Cartesian basis and consider a helix of radius r > 0 and pitch 27rb > 0 
parametrized by 

a(t) = r cos toi - 1 - r sint 02 + bte^ (SI.15) 


for t G M. Eliminating b between the relations 
Solving for b in the relation 

(f = arctan I - 


defining the lead angle and curvature of the helix yields 


and 


+ 52 


K, = 


COS^ (f 


(SI.16) 


(SI.17) 


If a plane strip approximated by a triangular lattice is wrapped helically around a cylinder, then at any point i, 
the curve in the j-direction is a helix. Further, letting be the lead angle of the helix through the point i in 

the j-direction, it is evident that 


+ I = + Y' 


(SLI 8 ) 


8 













Now substituting (SI. 18) into (SL17) yields 




K(2)(i) = 


.(3) 


(0 


cos^ (i) + |) 
r ’ 

COS^((/9^^^(i) + 


(SL19) 


Finally, substituting (SL19) into (SL8) yields 


Eb = ^ + cos"^ 0 

Reducing, as in the general case, results in 




'"<’>+1 


Eb = 


?>y/^kQA 

4^2 


(SL20) 


(SL21) 


Since = l/4 r^ for a helically wrapped strip, comparison of (SL21) with the Canham-Helfrich bending energy 
functional (SL5) confirms that 

(SL22) 


Notice that (SL21) is independent of (i) and thus of the orientation of the triangular lattice; consistent with the 
discussion of Vigliotti et al.^, this assures that the triangular lattice beha ves iso tropically in bending. Additionally, 
in the approximately unstretchable regime, where K is very close to zero, (SI. 12) is dominated by the term involving 
demonst ratin g that, for small k, we approximately minimize the integral o f over the surface. As reported 
by Sadowskyl ^ l ^ l ^ l, Mahadevan and Kellerl^ and Starostin and van der Heijdent^^^, this result is consistent with 
analogous results from continuum t heory. 

Finally, we can further compare (SI. 13) to the bending-energy density 


D{2H‘^ - (1 - v)K) 


(SI.23) 



and u = —1/3. 

8.2 Stretching energy 

Apart from energy due to bending, Mobius bands made from stretchable materials store elastic energy due to in¬ 
plane stretching. Using the purely two-dimensional theory of linearized elasticity and considering only neglecting 
contributions due to shear results in the area stretching-energy functional of the form (DesernoEl) 


-^stretch — 2^, 


L 


dA, 


(SI.24) 


where > 0 is the two-dimensional bulk (area) modulus and e ~ {J — l)/2 ~ {A/Aq — l)/2 is an approximation 
for the local in-plane strain, with J being the areal Jacobian determinant and A the current area of an area element 
with reference area Aq. (F^stretch carries dimens ions of energy.) In the present two-dimensional setting, ps carries 
dimensions of energy/length^. In view of (SI.24), the stretching-energy density is 


t^stretch — 2/ig6 


(SI.25) 


We now relate the energy of a discrete set of springs to the continuum expression (SI.24) with two primary 
objectives: (I) to determine the effective value of /i^ for an equilateral triangular lattice of linear springs, and (II) to 
verify the isotropy of such a lattice. The stretching energy Eg of a collection of linear springs is given by 




Ns 


Es = - rof 


(SL26) 


i=l 


where ki > 0 is the lattice modulus associated with stretching carrying the same dimensions as the two-dimensional 
bulk modulus Ng is the number of linear springs, r(i) is the current length of the i-th spring, and ro is the 
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equilibrium length of the springs. Notice that in the discrete setting, pure dilatation amounts to uniform stretching 
of the lattice. That is, each unit cell stretches (or contracts) into a larger (or smaller) equilateral triangle. The 
stretching energy of a unit cell in a periodic equilateral triangular lattice is 


Eve = ^[(?’(1) - rof + (r(2) - rof + (r(3) - rof], 


(SL27) 


where r(l), r(2), and r(3) are the lengths of the sides of the unit cells. Notice that the additional factor 1/2 in ( |SI.27 ) 
accounts for periodicity, namely that e ach spr ing is shared by neighboring unit cells. For pure dilatation (uniforrn 
stretching), r(l) = r(2) = r(3) = r and (SI.27) reduces to 


Eve = - rof. 


(SL28) 

Using the area Aq = v^rg/d of a unit cell and the magnitude of the pure dilatational strain e = (r — ro)/ro yields 

Es = TVuc^uc = dA. (SI.29) 

Js 


Importantly, (SI.29) is the lattice counterpart of (SI.24), and comparison of terms yields the connection 

V3ki 


(SI.30) 


between the effective continuum modulus /j^g and the lattice modulus ki. Notice that (SI.30|) is consistent with a 
previously establish connection between lattice and continuum parameters due to Seung and Nelson^^. 

The relationship ( SI.30|) is valid for periodic lattices. In the special case of one isolated unit cell without neighbors, 
the factor of 1/2 in ( |SI.27 ) vanishes, and it correspondingly follows that 


Msi 's/ski. 

Further, within the continuum theory, the following approximations hold 


tr (el) = 3e^J—1^— -1, 

^0 


(SL31) 


(SL32) 


where J is the areal Jacobian of the deformation and I is the two-dimensional identity tensor. In view of (SL32), 
the expression (SL29) for Eg can be written as 


Es = Nuc^uc 


kiN^crU A \ 

12 Uo ) 


(SI.33) 


8.3 Numerical results 



This occurs since, as bands approach a self-intersecting shape, the dimensionless continuum functional becomes 
orders of magnitude larger than the dimensionless lattice energy due to divergent values of and K near points of 
self-intersection (not shown). 

From Figure [SlT^ d), we observe that the dimensionless stretching energy due to dilatation is consistently of lower 
magnitude than that of the lattice model, both energy functionals exhibit the same trend. This indicates that for 
thin, stretchable materials capable of supporting shear strain, energy contributions due to in-plane shearing are not 
negligible when a strip of material is twisted into a Mobius band. As DesernoEl explains, the Canham-Helfrich 
theory neglects energy contributions due to shear and describes simple surfaces, like spheres, very well, but might 
not be suitable to model topologically complex objects such as Mobius bands. 
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Fig. SI.7: The dimensionless bending energy of the lattice model coincides with the bending energy of the Canham-Helfrich 
continuum functional using the theoretically determined splay modulus /x and saddle splay modulus /2: (a) Dimensionless 
bending energy of the lattice model (dashed lines) and that from the Canham-Helfrich functional (solid lines). Dimensionless 
energy contributions to the Canham-Helfrich energy functional from (b) and (c) K. (d) Dimensionless stretching energy 
due to dilatation (solid lines) is consistently of lower magnitude than that of the lattice model (dashed lines), but both exhibit 
the same trend. 


8.4 Numerical validation of the stretching energy 

For pure dilatation and assuming periodic boundaries, the total energy Eg of a collection of springs (SI.26|) and the 
energy of a unit cell are related through 

Es = A/'uc^uC: (SI.34) 


where A/'uc is the number of unit cells. In view of (SL34) and (SI.24), the effective stretching modulus can be 
extracted from the value of (SI.26) determined by simulation througE 




Ms 


2e^N\jcAo 


(SL35) 


By analogy to the variou s scal ings in troduced in the main text, using ki and the reference area A = NucAq to 
nondimensionalize (SI.34) and (SL35) yields the dimensionless counterparts 


4'., = 


E., 


E., 


Aki NuckAo 


and 


Ms 

Ms = ^ = 


E.. 


ki 2e^NucAoki 


(SI.36) 


(SI.37) 


In view of (SI.29), the following identifications hold: Es = a/3 dA, jig = a/3/2, and jlgi = a/3. To quantify 


the para meter range in which the plane lattice subject to in-plane pure dilatation behaves linearly elastically, we 
measure (SI.37) for different imposed e and different sizes of the latt ice (o r, equivalently, different A%c)- The lattice 
sizes range from 1 x 1 to 112 x 112 triangular unit cells (SI Figure SL8), and the strains range from e = —0.1 to 
e = 0.1. The energy is found to converge to the theoretically determined value, with little difference between results 
obtained for the 14 x 14 lattice and for the 112 x 112 lattice. 
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Fig. SL 8 : The lattice behaves linearly elastically, and the effective bulk modulus /Ig converges towards the theoretical predic¬ 
tion: (a) Dimensionless energy T 5 as a function of strain e along with the theoretical predictions T = \/3e^ and Ti = 2\/3e^. 
(b) Dimensionless effective bulk modulus /Ig as a function of strain e with the theoretical predictions /2g = \/3/2 and fig^ = ^/3. 
(b.i) Convergence of jUg for increasing number of segments n. 
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